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ABSTRACT 

Context. The idea of a unified description of supermassive and stellar black holes has been supported by the extension of the empirical 
radio/X-ray correlation from X-ray binaries to active galactic nuclei through the inclusion of a mass term. This has lead to the so-called 
fundamental plane of black hole activity in the black hole mass, radio and X-ray luminosity space. Two incarnations of this fundamental plane 
have so far been suggested using different underlying models and using two different samples of accreting black holes. 

Aims. We improve the parameter estimates of the fundamental plane and estimate the scatter of the sources around the plane in both samples. 
This is used to look for possible constraints on the proposed theoretical models. Furthermore, we search for selection effects due to the 
, inclusion of different classes of AGN or distance effects. 
Methods. We present revised samples for both studies together with a refined statistical analysis using measured errors of the observables. 
\Q | This method is used to compare the two samples and infer parameters for the fundamental plane in a homogeneous way. 
f*^) . Results. We show that strongly sub-Eddington objects in a state equivalent to the low/hard state of X-ray binaries follow the fundamental 
*** ; plane very tightly; the scatter is comparable to the measurement errors. However, we find that the estimated parameters depend strongly on 
, the assumptions made on the sources of scatter and the relative weight of the different AGN classes in the sample. Using only hard state 
I ■ objects, the fundamental plane is in agreement with the prediction of a simple uncooled synchrotron/jet model for the emitted radiation. 
O Inclusion of high-state objects increases the scatter and moves the correlation closer to a disk/jet model. This is qualitatively consistent with a 
picture where low-state objects are largely dominated by jet emission while high-state objects have a strong contribution from an accretion disk. 
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1. Introduction I2004I hereafter FKM). Similar unification efforts are also un- 

der way analysing and comparing the variability properties 



Active galactic nuclei (AGN) and black hole X-ray binaries of AGN and XRBs jyttlev et alJl2002t iMarkowitz et alJ l2C 

(XRBs) seem to have a similar central engine consisting of Kording & Falck el2004lAbramowicz et all2004l) . 
the central black hole, an accretion disk probably accompa- To establish connections between stellar and supermassive 

nied by a corona, and a relativistic iet Jshakura & Sunvaevl black holes we have to consider that black hole XRBs can be 

1223 Mirabel & Rodriguez! 1199% lAntonuccilll993h. Jet and found in distinct accretion states . In FKM we suggested that 



disk may form a symbiotic system jFalcke & Biermann||199fl a numbe r of AGN classes can be identified with correspond- 

|Falcke etalJh995|) which can be scaled over several orders ing XRB states and based on this proposed a power unification 

of magnitud e in mass and accretion rate llFalcke & Biermannl scheme for AGN and XRBs 

1996, 1999) suggesting that a single central engine can be used ™ . . , . vr)D , . , , 

1 ia ' " bb b b b The two most prevalent XRB states are the low/hard 

to describe very different types of black holes. tt /T u . . ... ■ , f . . . /uc . . 

J J ^ state (LH state) and the high/soft state (HS state, see e.g., 

While the general unification of stellar mass and super- iMcClintock & Remillard I2006I) . In the LH state the radio 

massive black holes picture has now been established for spe ctrum is alwa ys consistent with coming from a steady 

some time, it has recently been tested on a detailed empiri- jet (Fender 2001), which can sometimes be directly imaged 

cal level by correlations in the radio and X-ray band which <jStirlin^et^LTl200ll) . Once the source e nters the HS state, 

have l ead to the so-called f undamental plane of b lack hole ac- the radio emission seems to be quenched jFender et al. 1999; 

tivity jMerloni et alJl2003l hereafter MHDM. and lFalcke et alJ ICorbel etaD l2000l) . One possible scenario for the accretion 
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flow of a LH state object is that its inner part is opti- 
cally thin up to a transiti on radius , whe r e the flow turn s 
into a standard thin disk JEsin et al]ll997l IPoutanenl fT998). 
Usually, the X-ray emission of a LH state XRB is mod- 
eled using Compton ization (e.g. ISunvaev & Trumpeil 1 19791 
iThorne & Pricdll975h . however, some models suggest that the 
compact jet may contribute to the X-ray emission or even 
dominate it (IMarkoff et al.1 12001 [al IM arkoff & Nowakl l2004 
iHoman et al 120051 iMarkoff et all2005l) . As the diskfades. the 
system may become "jet-dominated" — meaning that the bulk 
of the energy o utput is in radiation and kinetic energy of the jet 

jFender et al 120031 FKM). 

Indeed, Corb el et al.ld2000l l2003) found a surprisingly tight 
correlation of the radio and X-ray fluxes of the black hole 
XRB GX 339-4 in its LH state which can be qualitatively 
and quantitati vely well understood in the c ontext of jet mod- 
els IMarkoff et al. 2003iL iGallo et alJ J2003I) showed that this 
correlation does not only hold for one source but seems to 
be universal for all LH state XRBs. It has also been observed 
that, once an object enters the high state, the radio emission is 

j ii — n 

quenched and it drops off the correlation ( Fender et al. 1999; 

lTananbau m£t^ljJl^2^JQ a ^ _£L a ^ 120031) . It has been sug- 
gested bv TMaccaroneetalTl2003^ a similar effect can be 
found in AGN. 

Radi o/X-ray correlations have a lso been found fo r AGN, 
e.g., bv lHardcastle & Worralll i 19991) and lCanosa etall lll999l) . 
The final breakthrough came when it was shown possible 
to combine these correlations to a fundamental plane in 
the radio/X-ray/black hole mass space for XRBs and AGN 
(MHDM, FKM). This fundamental plane gives a tight relation 
between the radio and X-ray fluxes and the black hole mass, 
which is valid for AGN as well as XRBs. Thus, the correlation 
proves the similarity of the central engines of these accreting 
black holes. 

However, there are at least two competing explanations 
for the fundamental plane. In the picture of FKM, the radio- 
through-X-ray emission for XRBs and the lowest luminosity 
AGN is attributed to synchrotron emission from a relativistic 
jet in the jet-dominated state (LH). As both components, the 
radio and the X-rays, originate from the same source - the jet 
- one can expect a tight correlation of both observables. We re- 
fer to this model as the 'jet only' model. One would expect the 
correlation to break down once a source leaves the radiatively 
inefficient accretion flow state and is no longer jet-dominated. 
Hence, the picture should not apply for high-state objects. 

On the other hand, MHDM suggested that the X-ray emis- 
sion originates from the accretion flow, while the radio emis- 
sion is still attributed to the relativistic jet. Both, the flow and 
the jet are presumed to be strongly coupled so that the radio and 
X-ray emission is correlated. Here we will mainly assume the 
accretion flow to be so me variant of an ad vection-dominated 
accretion flow (ADAF, Nar avan & Yil l 1994. se e also the con- 
vectio n dominated accretion flows e.g.. lOuataert & Gruzinovl 
2000) and refer to the model as ' ADAF/jet' model. The ADAF 
solution is only one possible acc retion flow model, e .g. on e 
other possibility is presented in lHaardt & Maraschil ll99ll) . 
Here we will use the ADAF/jet model only as the example for 
possible 'disk/jet' models. 



In the recent past, the statistics and slopes of the radio/X- 
ray/mass correlations have been u sed to argu e for and 
against the synchrotron/jet models llHeinzl Eoo4 MHDM). 
Hen ce, further cl arification is urgently needed. Additionally, 
iHeinz & Merlonil J2004I) have used the correlation to search 
for constrains of the relativistic beaming. However, as we will 
show here, all these analyses depend strongly on the statistics 
of the samples, the construction of the samples, and the as- 
sumptions on the scatter of the measurements. 

In this paper, we therefore investigate the problems of the 
parameter estimation of the fundamental plane of black hole ac- 
tivity. We will check the assumptions made by previous studies, 
and present an improved statistical analysis. We furthermore 
improve the samples presented by MHDM and FKM. With our 
refined parameter estimation method we analyze and compare 
both samples and investigate selection effects and the intrinsic 
scatter of the correlation. In this light, we discuss if the funda- 
mental plane can be used to constrain the underlying emission 
mechanism as previously suggested. We will use the intrinsic 
scatter to test which classes of AGN belong to the analog of the 
LH state XRBs. 

In Sect. 2 we discuss our method of parameter estimation, 
the improved samples, and discuss observing frequencies. In 
Sect. 3 we present our results and their implications and present 
our conclusions in Sect. 4. 

2. Parameter Estimation 

We are searching for the parameters of the fundamental plane 
for accreting black holes: 

logLx =£ R logL R +£ M logM + /? x , (1) 

where Lx is the X-ray luminosity in the observed band and Lr 
denotes radio luminosity at the observing frequency (vF v ), M 
is the black hole mass, the are the correlation coefficients, 
and bx denotes the constant offset. To simplify the notation we 
omit the units in the logarithms. Throughout this paper all lu- 
minosities are measured in erg/s, distances in pc and masses in 
solar masses. In the notation we follow FKM. To derive the pa- 
rameters as given in MHDM set £rx = 1/£ R and £ RM = £m/£r- 
The predicted values for the 'jet only' model are £r = 1.38 
and £ M = -0.81 (FKM), while the values for the 'ADAF/jet' 
model are £ R = 1 .64 and £ M = - 1 -3 (MHDM). 

2.1. The samples 

Here we compare the published correlations by MHDM and 
FKM. While both samples are used to extend the radio/X-ray 
correlation found in LH state XRBs, they differ in the selection 
of sources and which observing frequencies are used for the 
X-ray luminosities. 

MHDM sample 

The MHDM sample is a real radio/X-ray sample, i.e., it di- 
rectly uses the measured radio and X-ray fluxes. It contains 
XRBs and nearly all types of AGN except obviously beamed 
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sources like BL Lac objects. The sources were extracted from 
the literature under the condition that good mass estimates ex- 
ist. To obtain a representative sample the authors selected a 
similar number of bright active AGN and less active AGN. The 
sample contains low-luminosity AGN (LLAGN), LINERs (low 
ionization nuclear emission region), Seyfe rts (Type 1 and 2), 
FR Radio Gal axies JFanaroff ffi: Rileyl 19741) and radio loud and 
quiet quasars ( Kell ermann et alJ ll989) and the quiescence flux 
of Sgr A*. Thus, by using this sample one averages over nearly 
all types of AGN, whether they belong to the LH state or not. 

To avoid dealing with upper limits in the data we exclude 
all those limits from the MHDM sample. The overall result of 
the fit does not seem to change due to this, as we can repro- 
duce the best fit values of MHDM. This sample contains some 
XRBs that do not follo w the corre lation. Cyg X-l changes its 
state frequently JGallo et alJl2003l) and seems to stay always 
near the transition luminosity, thus, it will not trace the cor- 
relation well. GRS 1915+105 is a rather unique system that 
seems to stay in t he 'canonical' very high state most of its time 
(Rei g et alJl2003l) . We therefore exclude that source as well. It 
is still open whether LS 5039 is a black hole or a neutron star 
binary. Furthermo re, its radi o spectrum is peculiar for a LH 
state object jRibo et alJl2005l) . Thus, besides the original sam- 
ple we will also consid er the AGN sub sample of MHDM and 
add a subsample of the Gallo et al. ( 2003) XRB sample (see be- 
low). The quasar sample contains two very radio loud objects 
(3C273, PG 1226+023), while most quasars are radio quiet. To 
demonstrate the selection effects we exclude these two sources 
from the quasar sample when we consider subsamples of the 
MHDM sample. 

Our sample (KFC sample) 

The sample of FKM tries to include only objects in the LH 
state. FKM suggest to classify LLAGN, LINER, FR I Radio 
Galaxies and BL Lac objects as the analog classes of the LH 
state in XRBs. As FKM use a jet model as the basis of their 
suggested unification scheme, they try to compare observa- 
tions at frequencies that originate from synchrotron emission 
(see the discussion in Sect. 12.51 . They therefore extrapolate 
optical observations for FR I Radio Galaxies and BL Lac ob- 
jects to an equivalent X-ray flux (for details see FKM). The 
fluxes of the BL Lac sources have been deboosted with an 
average Doppler factor of 7 (FKM). In the current study we 
further increase th e number of LLAG N sources by including 
all sources of the iNa gar et alJ |2005) sample with Lr < 10 
erg/s for which we found X-ray fluxes in the literature. We 
will refer to this augmented sample as the 'KFC sample' 1 . 
The FKM LLAGN sample is ba sed on X-ray observa t ions o f 
the LLAGN sample studied by Terashima & Wilson (2003). 
Additionally we use fluxes from the followi ng surveys in or- 
der o f preference: The Chandra v3 pipeline JPtak & Gri ffiths 
2003 [www.x assist.org), the XMM serendipitous X-ray sur- 
vev JBarcons et alJl2002t). and the ROSAT HRI pointed cata- 
log JROSAT Scientific TearrfeOOOl) . Finally, NGC 4258 fluxes 

1 Despite other associations the reader may have with this abbrevi- 
ation, it is simply based on the present authorlist. 



were taken from lYoung & Wils^nl §004). In the KFC sample 
we on ly consider the non-Seyfert galaxies in the INagar et alJ 
(2005) sample; the Seyferts of this sample will be discussed 
separately as, even though they are of low luminosity, they still 
may belong to the supermassive analog of high state XRBs, 
because their black hole masse s are so low. For Sgr A*, we 
include the hard X-ray flare by Bagan offetalJ l200lh. as the 
flare may be due to jet emission (seee.g.. lMarkoffetalJ20 01b). 
Besides the flare we also show the result for the quiescent Sgr 
A* flux. 

XRB sample 

The sample of L H state XRBs is based on the sample of 
Gall o et alJ j2003). To avoid problems with state transitions 
we only include GX 334-9, V404 Cyg, 4U 1543-47, XTE 
1118+480, XTE J1550-564. For all sources we only con- 
sider the data if the source is in the LH state. We ex- 
cluded GRS 1915+105 and Cyg X-l as discussed above. For 
GX 334-9 JCorbel et alJ 120031) we used the updated X-ray 
fluxes from lNowak"etaI] J2005i) . 

2.2. Problems of Parameter Estimation 

The correlation between the radio and X-ray emission has been 
discussed for XRBs and AGN before and it has been shown 
by partial correlation analysis tha t the correlation is inde ed real 
(MHDM, for AGN only see lHardcastle & Worralll 19991) . For a 
discussion of the well constrained sources Sgr A*, NGC 4258, 
M81 and a XRB sample see Markoff (2005). Thus, we assume 
that the correlation exists and only check the parameter estima- 
tion process. 

To estimate parameters for measured data with err ors in 
all va riables, one has to use the merit function (see e.g.,[Press 
2002, MHDM). For measurements y; and Xij, which all have 
uncertainties, e.g., the y; can denote the X-ray luminosities 
while x\i denotes the radio luminosities and x^j the black hole 
masses, the merit function is defined as: 

„ 2 = y (yi-b- gj fljXi/ 

where the <x are the corresponding uncertainties, and the aj and 
b are the unknown parameters. The normal x 2 fits, which do 
not consider scatter in all variables, will yield asymmetric re- 
sults as for them only scatter in the "y"-axis is considered. The 
derivative of the merit function is non-linear in the aj and has 
to be solved with a numerical optimization routine. For the pa- 
rameter b one can still give an analytical solution: 

X>i(yi -bjXij) 
b mi „(a) = i- '—L (3) 

h w i 

with wr 1 = o-y + 2j (aj<x x . ; ) . Thus, we only have to solve a 
reasonably well behaved two dimensional function for which a 
standard numerical optimization routine can be used. 

The resulting parameters depend strongly on the assumed 
uncertainties cr in the data. MHDM assume that these uncer- 
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tainties cr are isotropic: 

o"l k = cr Lx = cr M . (4) 

Thus, they do not use measured uncertainties but set them 
isotropically to a value such that the reduced^- 2 is unity. This is 
a strong assumption and its effect has to be checked. As a first 
test we explore the effect of an anisotropy of the uncertain- 
ties in the mass estimation and the scatter in the luminosities, 
while still assuming that the uncertainties in the Lx-Lr plane 
are isotropic: 

o"L fl = cr Lx = 2ao-Q and cr M = 2(1 - a)cr . (5) 

With these assumptions of the uncertainties, we use the 
merit function to derive the parameters of the original MHDM 
sample. The best fit values are strongly depending on the 
isotropy parameter a as shown in Fig. An anisotropy pa- 
rameter of 0.5 corresponds to an isotropic distribution of the 
uncertainties. For this isotropic case we can reproduce the val- 
ues found by MHDM, which are also shown in the figure, a ~ 1 
corresponds to the case that the uncertainties in the luminosities 
dominate while for a « the uncertainties of the mass estima- 
tion are dominant. We observe that, for example, the parameter 
£r can take any value between 1 .4 and 3 for different a. 

In case that the uncertainties only deviate slightly from the 
isotropic case the slope of the parameters tells us how strong 
these errors propagate to the final fit values. Unfortunately the 
slope of the parameters around a « 0.5 is large. Thus, it is cru- 
cial to have a good estimate of the distribution of uncertainties. 
If we can improve the estimates of the uncertainties we can im- 
prove the validity of the parameter estimates. Any study based 
on a parameter estimate using isotropic uncertainties has to take 
the rather large additional uncertainties due to this assumption 
into account. 

Besides the demonstrated effect of the anisotropy of the un- 
certainties of the mass estimates and the luminosity estimates, 
a similar effect can be found if cr Lj) + ctl x We note that the best 
fit values do not depend on the absolute value of the combined 
cr but on the relative prominence of the different cr,. 

One problem of this parameter estimation scheme is that 
we can not deal with coupled uncertainties. We include several 
measurements of the XRBs GX 339-4 and V404 Cyg so the 
uncertainties of the mass and distance measurements are not 
independent for these datapoint. We will neglect this effect for 
simplicity, but it may influence the estimated parameters espe- 
cially for the small subsamples. 

2.3. Error Budget 

Both luminosities depend on the measured flux and the dis- 
tance, thus, the scatter in both quantities is coupled. To separate 
the equation into variables that have nearly independent errors, 
we separate the effect of the distance: 

Lr,x = F r ^r,xD 2 (6) 

where Frx denote the measured radio and X-ray fluxes and D 
denotes the distance. The E Rj x are conversion factors depend- 
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Anisotropy parameter alpha 

Fig. 1. The effect of anisotropic errors on the best fit parameters 
of the fundamental plane (eq.Q. The horizontal axis gives the 
anisotropy parameter a. A small value of a denotes that the un- 
certainties in the mass estimate dominate, while for large values 
the uncertainties of the luminosity are dominant, a = 0.5 corre- 
sponds to isotropic errors as used by MHDM. Also shown are 
the best fit values of MHDM and their uncertainties. Depending 
on the choice of the uncertainties one can obtain values of the 
parameters £r between 1.4 and 3 and £m between -0.8 and -3. 

ing on the observed band. We therefore find: 

log F x = £ R log F R + £ M log M + (2£ R - 2) log D + b x + log S, 

(7) 

where the mathematical conversion factors are combined into 
S. As they are just mathematical constants, they will not be 
discussed further. 

We have seen in the previous section that the assumption of 
isotropic uncertainties (as used by MHDM) will effect the best 
fit parameters in an unknown way. To improve this situation, we 
estimate the errors attributed to each variable. The fluxes Frx 
contain measurement errors and intrinsic scatter as discussed 
below. Besides this, the masses and distances are uncertain as 
well. However, as the correlation coefficient £r will be around 
1 .4 the effect of errors in the distance estimation are less severe 
for the correlation than for the individual luminosities as it ap- 
pears with a factor (2£r - 2). As the merit function requires 
Gaussian errors, we are always using symmetric errors in the 
log-log space, many of the errors below are indeed symmetric 
and the other parameters are only mildly asymmetric. 

- Mass estimate: Our XRBs LH state sample i s dom inated 
by GX 334-9 and V404 Cyg. For GX 339-4 lHvnes et alJ 
(2003) estimate a mass function of 5.8 + 0.5 M Q , which is 
therefore a lower limit for the mass of the black hole. We 
therefore assume a mass of 8 ± 2.0 M . For V404 Cyg we 
use a mass of 12 ± 2.5 M Q (Orosz (2003) gives a range 
of 10-13.4 Mq). Note that there are several data points for 
each object. Thus, their uncertainties are coupled, which 
we can not take into account. The mass es timates of the 
other XRBs are also taken from lbroszl J2003I) . For the AGN 
Merritt & Ferrarese boo lb give an absolute scatter of 0.34 
dex for M-cr relation. The mass estimate using the M-cr 
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relat ion is independent of the distance of the source (cf., 
Ferrarese & Ford 2005). This method is used for all galax- 
ies except our BL Lac objects. For masses estimated with 
the M-cr relation we will use the scatter of the correlation 
as a measure of the uncertainty for the mass estimate: 0.34 
dex. We use velocity dis persions from the Hypercat cata- 
log dPrugniel et al.lfl998l) . For the BL Lac obj ects we used 
indirec tly derived velocity dispersions from Woo & Urrvl 
(2002). Thus, these indirect measurements will have a 
higher uncertaintv. lBettoni et al 1 d200ll) give an uncertainty 
of these indirectly measured velocity dispersions <x a value 
of Scr = 18 km/s, which yields an additional uncertainty 
of « 0.3 dex for the mass estimate. Thus, we use an un- 
certainty of 0.46 dex for the mass estimate of the BL Lac 
objects. The indirect method to derive the velocity disper- 
sion depends slightly on the distance. However, as we only 
use this method for BL Lac objects and the distances are 
accurate compared to 0.46 dex uncertainty, we will ignore 
this effect. For Sgr A*, M 81, and NGC 4258 we use direct 
mass measurements (see FKM), the mass uncertainties are 
as low as 10% (< 0.05 dex). 

- Distance measurement: The dista nce o f GX 339-4 
is still under debate, (shahbaz et alJ d200 lh and 

I.Tonker & Nelemansl d2004h give a lowe r limit of 6 kpc, bu t 
the distance may be as high as 15 kpc (Hv nes et alJ l2004). 
We therefore adopt a distance o f 8 + 2 kpc. For V404 Cyg , 
we adopt the distance of 4 kpc djonker & Nelemanj2 004). 
We use an uncertainty of 1 kpc, as we can not account for 
asymmetric uncertainties and have the issue of coupled 
errors. For nearby AGN we us e upda ted distances from 

iMaoz et all ( 120051) : iTonrv et all d200ll) . If these are not 
available, w e use the distance estimates as given in the 

lTullvlill988h . The distance uncertainty is hard to access, as 
many different methods are used for which the uncertainty 
is sometimes not well known. We assume 40%, however, 
we checked that it does not change the result if one 
assumes less scatter. For AGN with distances derived 
from the Hubble law, we use an error estimate based 
on the peculiar velocities in the Hubble flow and the 
uncertainti es of the Hubble c onstant (we assume Ho = 72 
km/s/Mpc ISpergel et ai1l2003l 5% uncer tainty, fi A = 0.7 
and Qm = 0.3). iHawkins etal.1 d2003l) give a peculiar 
velocity of 506 km/s which corresponds to 6.7 Mpc. Thus, 
for most sources the distance uncertainty is mainly due to 
the uncertainty in the Hubble flow. For the MHDM sample, 
we use the distances as provided by MHDM and assume a 
constant uncertainty of 40% in the distance estimation to 
avoid that we overestimate the intrinsic scatter, see below. 
Here, we also checked that this assumption is not critical. 
For most sources MHDM derive the distance from the 
Hubble law, so that the uncertainty is mainly due to the 
Hubble constant. 

- Flux measurements: For all but the faintest objects are the 
fluxes very well constrained. Errors for radio and optical/X- 
ray fluxes are usually less than 10%. Systematic errors, e.g., 
due to the extrapolation of the different observed energy 
bands to our used X-ray band (0.5-10 keV), will also be of 
a similar magnitude. Thus, as these errors are small com- 



pared with the scatter due to the mass measurements, we 
do not introduce significant changes by assuming that these 
errors are isotropic. Due to this assumption, it is possible to 
unify these flux errors with the intrinsic errors of the source. 
- Intrinsic errors: Besides the measurement errors above, 
there are several sources of intrinsic scatter of often un- 
known magnitude: 

- Non-simultaneous observations of the AGN: The radio 
and X-ray observations are non-simultaneous, there is 
often more than a year between the different observa- 
tions. All accreting black holes (AGN and XRBs) are 
highly variable. Thus, already this effect can lead to a 
deviation by more than an order of magnitude. The ori- 
entation of this uncertainty is likely to be isotropic. 

- Beaming: In most models at least the radio emission 
is attributed to the relativistic jet and will thus be rel- 
ativistically beamed. In case that the X-rays originate 
from the disk/corona they will not be beamed and the 
deviations from the correlation will be enormous. For 
jet models, the X-ray emission may be beamed like 
the radio emission or have a different beamin g pat- 
ters (e.g., a velocity structu re in the ietlchiaberg e et all 
l2000l iTrussoni et ai1l2003l) . The asymmetry of this ef- 
fect depends on the exact model, so we can only assume 
isotropy. 

- Source peculiarities: The surrounding environment of 
the black hole will play a role on the exact emis- 
sion properties (e.g., there might be compact hotspots). 
There may also be an obscuring torus or other obsta- 
cles for the emission. This can result in strong X-ray 
absorption or the radio emission may also be absorbed. 
All models, however, only consider the nuclear emis- 
sion. 

- Spectral energy distribution (SED): Depending on the 
real emission model, it may be that we are not observ- 
ing the same emission type in the X-rays for the differ- 
ent objects. For jet models the effect of radiative cool- 
ing and the synchrotron cut-off have to be mentioned 
(see Sect. 12. 5\ . For disk models a similar effect may 
be due to the relative strength of the disk component, 
the jet component as observed in the resolved X-ray 
jets, and the Comptonization component. The X-rays in 
XRBs may not originate from the same process as those 
in AGN due to the mass scaling by whatever theory is 
used. 

The total intrinsic scatter will be derived from the scatter of 
the correlation. We will see below that the intrinsic scatter 
is surprisingly small considering this long list of possible 
errors. 

For the two last sources of scatter, the flux measurements and 
the intrinsic errors, we do not have an exact knowledge of their 
magnitude and their asymmetry. We will therefore assume that 
they are isotropic in the log Fx - log plane and parameterize 
their combined magnitude as cr In ,. This parameter <j\ nt will be 
chosen, such that the reduced merit function is unity. With this 
choice we assume the error distribution just described; this will 
therefore affect the final fit values. If one distributes the excess 
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variance in a different manner one will find slightly different 
best fit parameters. However, as we know several effects that 
introduce uncertainties in the radio and X-ray fluxes and the 
the excess variance in these variables is surprisingly low, it is 
sensible to include the excess variance only in the radio and 
X-ray fluxes. 

2.4. The parameter cr\ nt 



Lg F e 



In Sect. |221 we present our assumptions on the uncertainties of 
the measured variables. In case that these assumptions are ex- 
act, e.g., the magnitude of the uncertainties and that they are 
Gaussian distributed, the fitted parameter cr In , will describe the 
real intrinsic scatter of the sources. This value could then be 
used to constrain the different contributions like the effect of 
beaming. However, if we overestimate the uncertainties in the 
flux, distance, and mass estimation the derived cr Inf will be too 
small. Similarly, an underestimation will lead to an overestima- 
tion of the intrinsic scatter. 

Our method can not treat coupled uncertainties correctly. 
For XRBs we include several data points for each source, but 
there is only one mass and distance estimates for that source. 
Thus, as we have to assume that our uncertainties are inde- 
pendent, we overestimate the measurement errors. Thus, the 
inferred cr In , for XRBs alone is zero. This means that the devi- 
ations of the data points from the optimal correlation is within 
the measurement uncertainties. On the other hand, for the AGN 
samples the uncertainties are independent. Every source is in- 
cluded only once in the sample. Here, cr In , should be a good 
measure for the intrinsic scatter. 

2.5. Origin of X-ray emission in the jet model 

The spectrum of a relativistic jet can be directly observed 
in BL Lac objects, as relativistic boosting increases the rel- 
ative prominence of the J et component compared to the disk 
(Blandford & Reeslfl978l) . An idealized spectrum of a jet, i.e., 
the "Camel's back", is shown in Fig.[2]in flux (F f ) representa- 
tion. Such a jet component exists at least in every AGN with a 
detectable jet, most likely in all AGN. The relative prominence 
of this component in respect to disk and corona emission will 
vary. Furthermore, the exact shape of the SED depends also on 
the inclination angle, the Lorentz factor of the jet and peculiar- 
ities of the source. 

The spectrum of a conical jet is flat (in F v representation) 
due to optically thick synchrotron emission up to the turnover 
frequency e t . This is followed by the optically thin power law 
component with a typical energy index a (F e ~ e~ a ) of around 
a = 0.5. This emission comes from the innermost region of the 
radiating jet. The power law continues up to the energy where 
radiative cooling plays a role (e coo i); here the power law may 
steepen to a as 1. Finally the synchrotron emission cuts off 
at (e cu i), due to the acceleration mechanism creating the radi- 
ating particles. At photon energies above the synchrotron cut- 
off, synchrotron-self Compton emission and external Compton 
emission are visible. As the jet is relativistically boosted, the 
total observed power in these two humps depend on the jet 



Synchrotron-Self Compton/ 
External Compton 




e cool E cut 



Lg photon energy 



Fig. 2. Sketch of the SED of a relativistic jet as for example 
observed in BL Lac objects. At lower photon energies syn- 
chrotron emission is dominant while the high energy spectrum 
is usually explained by inverse Compton processes. See text for 
details. 



Lorentz factor and the inclination of the source. The relative 
prominence of the synchrotron and the inverse-Compton emis- 
sion depends on these two parameters as well. 

The photon energy where the synchrotron emission cuts off, 
seems to depend on the total power of the accreting system, see 
e.g., t he blazar sequence JFossati et alJ ll998: Ghiselli ni et all 
2002). For low power systems, like high-peaked BL Lac ob- 
jects (see the previous references), the cut-off energy can be 
above or around the standard X-ray band of 0.5-10 keV, while 
for strongly accreting systems, e.g., flat spectrum radio quasars, 
this cut-off can be as low as 0. 1 eV. Thus, we can expect that 
only the X-rays of low power systems actually originate from 
synchrotron emission. 

Radiat ive co oling in the X-rays, as for example discussed 
by|Heinz (2004), may play a role in some sources of low to 
intermediate accretion rates. Its main effect will be that the 
measured X-ray flux will be reduced to what one would expect 
from a simple uncooled jet. However, the X-ray reduction due 
to cooling is far less severe than the synchrotron cut-off. The 
effect of cooling will increase the observed intrinsic scatter as 
we are exploring only the simplest jet model: the uncooled jet. 
For most AGN besides the LLAGN in the MHDM sample this 
synchrotron cut-off will be below the X-ray band. The X-ray 
emission in these sources will not be due to synchrotron emis- 
sion. It will either originate from inverse-Compton processes in 
the jet or from the disk/corona. 

The cut-off energy (ecut) of LLAGN sources in the KFC 
sample will for most sources be above the ROS AT and Chandra 
bands. The radio luminosities of the LLAGN are in the ra nge of 
2Q36-38 er g/ s if one extrapolates Fig. 7 of lFossatiet 
to these energies, one can expect a cut-off above 10 19 ~ 20 Hz or 
40-400 keV. For BL Lac objects and FR-I RGs, the KFC sam- 
ple extrapolates optical fluxes to an equivalent X-ray flux. This 
treatment avoids the effect of the synchrotron cut-off also for 
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these sources. The 'jet only' model should therefore be appli- 
cable to the KFC sample. 

2.6. Origin of the X-ray emission in the MHDM sample 

The X-ray emission of many sources in the MHDM sam- 
ple is very likely not due to synchrotron emission, due to 
the synchrotron cut-off. For example, for radio loud quasars, 
the synchrotron cut-off is far below the X-ray regime (see 
e.g.. lTavecchio et afll2002l) . Only the X-ray emission from the 
LLAGN may originate from synchrotron emission. Thus, the 
'jet-only' model is not applicable to the MHDM sample, and 
can therefore not be constrained using the MHDM sample. 

The AGN sample of MHDM contains radio loud objects 
(e.g., Cyg A, 3C273, etc). The X-ray spectral index of radio 
loud Quasars (RLQ) and those of radio quiet Quasars (RQQs) 
seems to be different: Fr^q « 1 .63 + 0.02 compared to Frqq a 
1.89 + 0.11. Furthermore, RLQ show very weak or no reflec- 
tion components and the strength of the soft excess seems to be 
anti-correlated with the radio-loudness. These effects are usu- 
ally explained by a relativistic jet component in the R LQ case 
(see e.g.. lReeves & TurneifcOOotlPiconcelli et al.l2 005). 

On the other hand, there are several detected "resolved" X- 
ray jets in AGN. For th e source 3C273, whic h is included in 
the MH DM sample, see Marshall et alJ J200ll) . iMarshall etail 
J2005h and ISambruna et alJ d2004 find in more than 50% 
of the observed radio loud sources resolved X-ray jets. The 
authors suggest that all bright radio jets may have X-ray 
counterparts. The X-ray emission from the AGN jets is usu- 
ally explained by non-thermal processes (synchrotron, inverse- 
Compton, synchrotron-self-Compton), which process domi- 
nate s seems to vary from source to source. For a discussion 
see lHarris & Krawczvnskil J2002). As the jet is visible, it most 
likely contributes at least at some level to the core X-ray flux. 
Furthermore, many of the MHDM fluxes are derived from 
ASCA or GINGA data. In that case, the Chandra resolved X- 
ray jet will be observed as a point source. Note that the X- 
ray jets in AGN are often dominated by emission from knots, 
unli ke what i s expec ted for LH state XRBs as modeled by 
iMarkoff et alJ d2001ah . From the statistical studies and the di- 
rect observation of X-ray jets, we conclude that at least some 
radio loud AGN in the MHDM sample have an X-ray compo- 
nent originating from the jet. 

Besides jet components, other features can contaminate the 
X-ray fluxes, especially for non-Chandra data. The flux may 
consist of several components originating from different phys- 
ical processes, including the disk described by any model, the 
corona, the reflection component, warm gas and the jet (at least 
in radio loud objects). 

2.7. Uncertainties of the estimated parameters 

Given the data points and their estimated uncertainties we can 
derive the optimal fit parameters from eq. |2] The merit func- 
tion x 2 can be used similarly to the usual x 1 to estimate the 
confidence region of the parameters. The 1 <x confidence re- 
gion should be given by A£ 2 « 2.3 as we have 2 degrees of 



freedom in our model besides the offset bx- As problems may 
arise due to the use of the nonlinear merit function or the un- 
known distribution function of the errors, we checked that this 
confidence region is in agreement with the confidence region 
derived by a Monte Carlo simulation and the Bootstrap method 
(see below). 

In the Monte-Carlo simulation, the errors of the parameters 
are estimated by creating a large number (5000) of artificial 
datasets that have similar statistical properties compared to the 
measured dataset. For each of the artificial datasets, we esti- 
mate the best-fit parameters using the same method as for the 
original dataset. From these fitted parameters, we derive the 
confidence region and the A£ 2 corresponding to the 1, 2 and 
3 cr confidence regions. To create the artificial data, we con- 
sider each population of sources of our measured dataset indi- 
vidually and measure their scatter and position in the log Lr - 
log Lx plane. In the artificial dataset, we distribute the sources 
uniformly over their radio luminosities. The artificial dataset 
will therefore also contain simulated objects of all considered 
types. We assume that the scatter compared to the real correla- 
tion is Gaussian. As seen in Fig.[3]this seems to be roughly the 
case. We also introduce flux limits in the simulation to access 
the effects of the distance selection effect. 

The bootstrap method functions as follows. From a set of 
N measured sources, we draw N at random with replacement, 
thus, creating an artificial dataset. This dataset contains some 
of the sources more than once, while others are omitted. The 
parameters of the fundamental plane will be estimated for this 
sample in the same way as for the original dataset. These simu- 
lated parameters should be distributed around the original best 
fit values as the measured pa rameters are distributed around the 
real parameters ( Press 2002). The benefit of this method is that 
it does not require prior knowledge of the distribution function 
from which the original dataset was drawn. 

2.8. Different parameter estimators 

Up to now we have only discussed the parameter estimation 
using the Merrit function (eq. This method minimizes the 
average distance of the data-points from the plane weighted 
with the measured uncertainties. Other possible methods used 
in astronomy include the ordinary least sq uares estimation o r 
Maximum Likelihood methods (see e.g., D'Agostini 2005). 
Even though Maximum Likelihood estimators find the most 
probable parameters of a model, the method is often biased to- 
wards lower fit parameters. To find the optimal fitting method 
for our problem, we compare the different fitting methods using 
a Monte Carlo simulation. We create several artificial samples 
with our Monte Carlo simulation and compare the results of 
the different estimators with the parameters used to create the 
sample. We set the intrinsic scatter of our artificial XRB and 
LLAGN sample to 0.2 and for BL Lac and FR-I RGs to 0.8, 
which is roughly double of what is found in our sample. The 
resulting intrinsic scatter of the simulated sample is cr,„, = 0.65. 
For each parameter estimator consider two cases: First, we only 
use the average intrinsic scatter (0.65) to estimate the parame- 
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Table 1. Performance of the different parameter estimators. The given uncertainties is the standard deviation of the fit parameters 
obtained from different artificial datasets. 



Estimator 










X 


Assumed Parameters 


1.40 


-0.85 


-4 


.9 


Merrit Function 


1.42 ±0.13 


-0.87 ± 


0.16 


-5.4 


±3.6 


Merrit Function with exact knowledge of scatter 


1.40 ±0.08 


-0.85 ± 


0.09 


-4.8 


±2.3 


Merrit Function with isotropic uncertainties 


1.54 ±0.16 


-1.03 ± 


0.19 


-8.5 


±4.4 


Maximum Likelihood 


1.10 ±0.08 


-0.47 ± 


0.09 


3.5: 


t2.1 


Maximum Likelihood with exact knowledge of scatter 


1.25 ±0.06 


-0.68 ± 


0.07 


-0.51 


± 1.68 


Ordinary least squares 


1.09 ±0.08 


-0.46 ± 


0.09 


4.0: 


t2.3 



35 
30 
25 
20 
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5 




LH sample 
Gaussian 



Li 



— 



-4 -3-2-101234 
Deviation from fit 

Fig. 3. Histogram of the scatter in the full KFC sample. The 
scatter can be well approximated by a Gaussian with the 
standard-deviation cr = 0.50. 

ters, and second, we use the exact probability distribution used 
to create each individual data-point. 

For each estimator we simulated 100 different datasets and 
give the average estimated parameters and the standard devia- 
tion in table !2.8l While the Merrit function seems to be a fairly 
robust method, the maximum likelihood estimator is biased to- 
wards smaller fit values. Thus, we will estimate our parameters 
with the method using the Merrit function described in section 

El 

3. Results 

Our fitting algorithm is only robust if our sources are normally 
distributed around the fundamental plane. As a first test we 
show a histogram of the scatter of the KFC sample around the 
best fit to the fundamental plane in Fig. [3] The deviations are 
roughly normally distributed with cr x 0.5. Thus, the developed 
analysis method can be utilized. 

We fitted the KFC sample with the fundamental plane de- 
scribed in eq.0 The parameter estimation method is described 
in Sect. |2 We find as best fit values 



= 1.41 



+0.14 
-0.12 



-0.87 



+0.15 
0.17 



-5.01 



3.9 



(8) 



The confidence region of the two relevant parameters and the 
fit is shown in the top panel of Fig.0] The given uncertainties 



are derived from the A£ 2 map. The intrinsic scatter in the fun- 
damental plane is for this sample <j m = 0.38+0.06. The uncer- 
tainty of (Tim has been derived by bootstrapping the measured 
sample. 

The results of the KFC sample and the MHDM sample and 
their subsamples is presented in Table [2] The errors given in 
this table are derived using the Bootstrap method. As the er- 
rorbars for the KFC sample are nearly identical for both meth- 
ods, we present only the bootstrapped errorbars for the different 
subsamples. This error estimation is less dependent on the as- 
sumption of normal distributed scatter. As the uncertainties are 
not strongly anisotropic, we give a symmetric uncertainty to 
avoid problems with anisotropic errors. We note that the differ- 
ent fit values do not coincide within the errors for both samples 
due to the different source populations included in the sample. 

The subsample of the KFC sample containing no FR-I ra- 
dio galaxies has a larger radio coefficient (^r = 1.64 + 0.13) 
and a smaller mass coefficient (^m = -1.08 ± 0.14) than the 
full sample. On the other hand, the subsample without BL Lac 
objects deviates in the other direction (£r = 1.25 ± 0.1 and 
£m = -0.74 ± 0.12). This difference arises, because FR-I ra- 
dio galaxies are fainter in the optical and X-ra ys than the de- 
boosted BL Lac objects ( Chiaberg e et al. 2000).' The difference 
in the optical/X-ray luminosity can be explained by a velocity 
structure in the jet. As we include a similar number of BL Lac 
objects and FR-I radio galaxies in the full sample, these effects 
will partly average out. However, the final value for the param- 
eters will depend on the total weight each AGN class has in the 
total sample. 

The KFC subsample containing only XRBS, Sgr A* and 
LLAGN has similar best fit values as the subsample contain- 
ing no FR-I RGs. Its correlation coefficient, £r = 1.59 ± 0.21, 
is larger than the one found for the full sample. This is partly 
due to the fact that the intrinsic scatter of this subsample is 
extremely low: cr ; „, = 0.12. Thus, the errors of the mass esti- 
mation dominate (0.34 dex) the overall error budget and one 
gets larger fit values as shown in Fig. ^ If one adds Seyferts 
and Transition objects, which may correspond to the high state, 
the correlation coefficient gets even larger. However, all KFC 
subsamples seem to be roughly in agreement with ^r ss 1 .4 
and £m ~ -0.8. This radio coefficient is in agreement with 
the radio/X-ray cor relation for GX 339-4, which has ^ * 1.4 
JCorbelet all2003l) . 
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Table 2. Best fit results for the KFC sample and the original and edited MHDM sample. Besides the full sample we also give the 
parameters for the subsamples containing only a limited set of AGN classes. Note that LINER sources classified by MHDM are 
not limited to LLAGN. The LH XRB sample is defined in Sect. 12. II The column N denotes the number of sources in the sample. 





ft 






bx 


(Tint 


tm)gr= 1.4 


(7to£r= 1.4 


N 


KFC sample 


Full Sample 

(XRB, Sgr A*, LLAGN, FR-I, BL Lac) 


1.41 ±0.11 


-0.87 ±0.14 


-5.01 ±3.20 


0.38 


-0.86 ± 0.02 


0.38 


77 


Full Sample + Seyferts & Transition obj. 


1.48 ± 0.13 


-0.95 ±0.16 


-6.89 ± 3.83 


0.44 


-0.86 ± 0.02 


0.44 


100 


XRB, Sgr A*, LLAGN, FR-I 


1.25 ±0.10 


-0.74 ± 


0.12 


-0.46 ± 2.93 


0.28 


-0.91 ± 0.02 


0.30 


58 


XRB, Sgr A*, LLAGN, BL Lac, 


1.64 ±0.13 


-1.08 ± 


0.14 


-11.67 ±3.59 


0.18 


-0.81 ±0.02 


0.23 


62 


XRB, Sgr A*, LLAGN 


1.59 ±0.21 


-1.02 ± 


0.21 


-10.15 ±6.17 


0.12 


-0.84 ± 0.02 


0.15 


43 


XRB, Sgr A*, LLAGN, Seyfert & Transition 


1.86 ±0.35 


-1.33 ± 


0.36 


-17.98 ±9.92 


0.35 


-0.85 ± 0.02 


0.39 


66 


Sgr A*, LLAGN, FR-I, BL Lac 


1.70 ± 1.17 


-2.17 ± 


4.09 


-5.21 ± 12.08 


0.45 


-1.15 ±0.38 


0.46 


52 


Full Sample with quiescent Sgr A* 


1.52 ±0.14 


-1.00 ± 


0.18 


-8.15 ±4.11 


0.39 


-0.85 ± 0.02 


0.40 


77 


XRB, quiescent Sgr A*, LLAGN 


1.97 ±0.11 


-1.41 ± 


0.10 


-21.05 ± 3.05 


0.00 


-0.86 ± 0.02 


0.23 


43 


Original MHDM sample 


1.45 ±0.17 


-0.99 ± 


0.22 


-5.98 ± 5.02 


0.72 


-0.93 ± 0.03 


0.73 


116 


MHDM Sample with LH XRBs 


Full Sample 

(XRB, Sgr A*, LINER, Quasar, Seyfert) 


1.74 ±0.20 


-1.35 ± 


0.27 


-14.23 ±5.75 


0.65 


-0.92 ± 0.03 


0.68 


103 


XRB, Sgr A*, LINER, Seyfert 


1.55 ±0.19 


-1.15 ± 


0.24 


-8.93 ± 5.37 


0.64 


-0.96 ± 0.03 


0.64 


92 


XRB, Sgr A*, LINER, Quasar 


2.12 ±0.31 


-1.75 ± 


0.38 


-25.20 ±8.81 


0.51 


-0.90 ± 0.04 


0.62 


57 


XRB, Sgr A*, Quasar, Seyfert 


1.79 ±0.28 


-1.39 ± 


0.37 


-15.68 ±7.87 


0.63 


-0.88 ± 0.03 


0.65 


82 


XRB, Sgr A*, LINER 


1.59 ±0.35 


-1.19 ± 


0.40 


-10.11 ± 10.13 


0.50 


-1.54 ±0.21 


0.71 


46 


XRB, Sgr A*, Quasar 


2.03 ± 0.28 


-1.55 ± 


0.36 


-22.76 ± 7.86 


0.14 


-0.74 ± 0.04 


0.32 


37 


Sgr A*, LINER, Quasar, Seyfert 


1.65 ±0.26 


-1.72 ± 


0.50 


-8.02 ± 8.06 


0.74 


-1.41 ±0.22 


0.76 


78 



If one uses the quie scent flux of Sgr A* instead of the flare 
by Baganoff et al I fcOOlh in the KFC sample the best fit values 
for^R increase. The subsample containing only XRBs, the qui- 
escent Sgr A* and LLAGN yields £ R = 1.97 and £ M » -1.41. 
Interestingly, the latter subsample has an intrinsic scatter of 
zero, i.e., the scatter of the correlation is in agreement with the 
assumed measurement errors. However, the inferred fit for £r 
is no longer in agreement with the value found for XRBs only 
(£r = 1 .4 ). The quiescent X-ra y emission from Sgr A* is an 
extended (Bagano ff et aljE003l) . while this is not the case for 
cores of the XRBs or AGN. Thus, it is not surprising that the 
fit values change. 

For the previous fits we have used the same intrinsic scat- 
ter (Ti„, for all our objects. However, we have seen that XRBs 
and LLAGN can be fitted with significantly less excess scat- 
ter than the full sample. Thus, we can fix the intrinsic scatter 
of XRBs and LLAGN to 0.1 dex and fit only the scatter in the 
BL Lac and FR-I RGs. Now, we find ^r = 1.61 ± 0.11 and 
£m = -108 which is similar to the result found for XRBs, 
Sgr A* and LLAGN only. The intrinsic scatter found for FR- 
I RGs and BL Lac objects is 0.55 dex. The statistical weight 
for these two classes is therefore significantly less than that of 
the LLAGN and XRBs. It is therefore not surprising that we 
find similar fit values: the fitting method puts only very a low 



weight on the additional sources. To observe selection effects 
it is therefore sensible to use a constant intrinsic scatter for all 
sources. 

For the full edited MHDM sample, we find £ R * 1.74+°-^ 
and £m ~ - 1 -35^030. As f° r tne KFC sample, we also find 
a selection effect for the MHDM sample. The MHDM sub- 
sample containing no Seyfert objects yields as best fit param- 
eters £r « 2.12 and £m ~ -1.75, while the one containing no 
Quasars gives ^r » 1.55 and £m ~ -1.15. This effect can be ex- 
plained by the fact that the quasars are more radio quiet than the 
other AGN in the sample. We note that the formal fits of £r for 
most subsamp les tend t o be s ignificantly higher than the value 
found bv lCorbel et ail d2003l) of ^r = 1.4 for XRBs. Overall, 
we find the choice of the sample strongly influences the final fit 
value in both considered samples (MHDM and KFC). This is 
not necessarily worrisome, as one does expect somewhat dif- 
ferent results for different black hole states. 

The correlation index £r for XRBs only is better con- 
strained than in the case we are considering here, as one does 
not have to consider the mass scaling in that case. For XRBs 
only one finds & * 1.4 dCorbel et alJl2003UGallo et alJl2003h . 
If we use this prior knowledge to fix the correlation index £r 
to 1 .4 and only fit the mass scaling parameter £m to our full 
KFC sample, we find £ M = -0.86 + 0.02 (see table|2). This fit 



10 



Elmar Kording et al.: Refining the fundamental plane 



-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 
-1.6 
-1.1 

-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 
-1.6 
-1.1 

-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 
-1.6 
-U 



Jet-Model 
Jet/ADAF 



1.4 1.6 1.! 

Radio Exponent 



Jet-Model 
Jet/ADAF 



1.4 1.6 1.! 

Radio Exponent 



Jet-Model 
Jet/ADAF 



1.2 1.4 1.6 1.1 

Radio Exponent 




Jl 50 

i 

X 

. I 45 



Full Sample + 
Best Fit 



28 30 32 34 36 38 40 42 
Log L_R [erg/s] 





Only LLAGN+XRB + 




Best Fit 













28 30 32 34 36 38 40 42 
Log L_R [erg/s] 





Edited MHM sample +**' 
Best Fit ;*-V— 




+ +,'' * 













28 30 32 34 36 38 40 42 
Log L_R [erg/s] 



Fig. 4. On the left side we show the \ 2 map with one (solid line), two (dashed) and three sigma (dotted) levels together with the 
predictions of the 'jet only' model and the ADAF/jet model. The right side shows the best fit and the AGN and XRB sample. The 
three rows show the different samples, from Top to bottom: KFC sample, Only XRBS and LLAGN of the KFC sample, and the 
edited MHDM sample. 



is shown in Fig. [5] The fit is not significantly worse than the 
original fit in both parameters, as the intrinsic scatter increases 
by less than 0.01 (cr m , - 0.38). Interestingly, the best fit value 
for the subsample of LLAGN, that had a higher correlation 
(£ R = 1.59) yields with the prior knowledge f M = -0.84 + 0.02, 
in agreement with the value found for the full sample. For the 
MHDM sample, we find £ M = -0.92 + 0.04. As the correlation 
index £r is well constrained for XRBs, fixing £r = 1 .4 will 
likely yield the best estimates for the parameters of the funda- 
mental plane. 

Besides checking the effect of anisotropic errors we also 
explore the effect that the uncertainties are not normally dis- 
tributed but exponentially. In that case, one does not mini- 
mize the square of the deviations divided by th e uncertainty 
but the absolute value. This "robust" estimation JPressl Eo02) 
yield £ R = 1.39 and f M = -0.84. For the edited MHDM sam- 
ple we find £ R = 1.63 and £ M = -1.14. Thus, the change of the 



distribution of the uncertainties does not significantly change 
the best fits. 

In table |2] we also show the intrinsic scatter measured for 
the different datasets. If we drop our assumption that the intrin- 
sic scatter is the same for AGN and XRBs, e.g. if cr ml is smaller 
for XRBs than for AGN, then the correlation coefficient £r will 
rise slightly. 

3. 1 . The distance selection effect 

Our sample of AGN and XRBs is a random sample of well 
studied objects and not a complete distance limited sample. 
Thus, one can fear that the fundamental plane is a spurious 
correlation created by distances selection effects. It has been 
shown with partial correlation analysis by MHDM that the fun- 
damental plane is real and not spurious. Th is problem has been 
discussed in detail bv lMerloni et alJ (|2006), where the authors 
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Fig. 5. Best fit for the full KFC sample if one fixes ^r to 1 .4, 
which is the value found for XRBs. Only the mass scaling pa- 
rameter is fitted. One finds: £ M = -0.86 + 0.02 and b x = -4.9. 
The intrinsic scatter increases by less than 0.01 compared to 
the fit in both parameters (en,,, = 0.38). 



fluxes and do not include upper limits which could affect the 
significance. 

On the other hand, if we fix = 1 -4 the fundamental plane 
suggests that -M- and M are correlated. The Kendall r for this 

correlation is r = 0.6, which is significant as for uncorrected 
data t is normally distributed around with a standard devia- 
tion of 0.006. Again we can check whether this is due to the 
different distances in the sample: The partial correlation coeffi- 
cient is r Lx/L iA MD = 0.84, which is again significant (7 <x). 

Even if the observational flux limits can not create a spu- 
rious fundamental plane, these limits might still bias the esti- 
mated parameters. This problem can be tested with our Monte 
Carlo simulation. We start with our fundamental plane and ob- 
serve the parameter changes due to increasing flux limits. As 
the flux limits may reduce the observed scatter around the fun- 
damental plane, we increase the intrinsic scatter by a factor of 
2. The results are summarized in table [5] As long as the flux 
limits stay in a reasonable range, the changes to the parameters 
are below 0.1, i.e., they are of the order of the uncertainties of 
the fits. 



show further statistical and observational evidence that the fun- 
damental plane is indeed real. Here, we will just present some 
further tests for our sample. 

The fundamental plan e is an exten sion of radio/X-ray 
correlations in XRBs (Gall o et al.l 120031) . This correlation is 
highly significant and is not effected by distance effects as 
we can trace individual obj ects on the correlation. I n AGN, 
lHardcastle & Worralll Jl999) and ICanosa eTaD dl999l) showed 
that a radio/X-ray correlation exists and is also not created as 
an artifact of plotting "distance against distance" in a flux lim- 
ited sample. If one combines thes e two correlations wi th a mass 
term (which is needed, see Falcke & Biermann 1996j), one ar- 
rives naturally at the fundamental plane. 

One method to check if a correlation between two observ- 
ables is only created by a third variable, is the partial correlation 
coefficient. In our case for radio and X-ray luminosity is the pa- 
rameter that may create a spurious correlation the distance. The 
correlation coefficient is defined as: 



' rx,d — , » \y i 

^(l.-^xi.-^) 

where r a b is the normal Pearson R for a and b. This correlation 
coefficient should remove the effect of the different distances. 

If we fix the mass coefficient £m = -0.85, we find r rx j = 
0.91. We have written a small Monte Carlo simulation, to test if 
this coefficient is significant. We assume a uniform distribution 
for the distances and the fluxes in log-space, and assume that 
the fluxes are uncorrected. We created 10 6 artificial datasets, 
and found that not a single dataset had a partial correlation co- 
efficient as large as 0.8, even though the average normal cor- 
relation coefficient for the luminosities was 0.99. The standard 
deviation of the partial correlation coefficient is 0.11 and the 
mean, 0. Thus, the radio/X-ray correlation with a fixed mass 
coefficient is significant at the 8 cr level. However, in our sam- 
ple, we only include objects with measured radio and X-ray 



3.2. Comparison of the KFC sample and the MHDM 
sample 

The correlation coefficient £r in the MHDM sample and its 
subsamples seems larger than those usually found in the KFC 
sample and a similar effect can be found for £m- The fits of 
the KFC sample and its subsamples are in agreement with the 
value found for LH state XRBs (£r = 1 .4), while the deviations 
are larger for the MHDM sample. This may be seen as a hint 
that the MHDM sample is not simply a continuation of the LH 
state correlation for XRBs, but it may contain other effects like 
a different source of emission. 

The main statistical difference between the MHDM sample 
and the KFC sample is the smaller intrinsic scatter of the latter: 
(Ji n , = 0.39 compared to <j m , = 0.65. It is hard to assess the 
uncertainties of this value due to selection effects as the under- 
lying distribution is unknown. Bootstrapping yields an error for 
both values around 0.06. 

The discrepancy is partly due to the fact that the KFC sam- 
ple has less AGN compared to XRBs than the MHDM sam- 
ple. Ho wever, even if one a dds the Seyferts and transition ob- 
jects of Naga r et all d2005h to create a sample of similar size 
than the MHDM sample, then the intrinsic scatter is still less 
((Ti nt = 0.46) than for the MHDM case. The same is true for the 
subsamples of similar size. 

The most homogeneous subsample is the sample contain- 
ing only LLAGN, LH state XRBs and Sgr A*. Here we find 
cr m , — 0.1 1. This low value is not only due to the fact that we 
overestimated the errors of the XRBs, as the numerical value 
is below 0.1 for the LLAGN sample without XRBs as well. 
Thus, the correlation is extremely tight for the lowest luminos- 
ity objects. If we extend this sample to slightly higher accretion 
rates, i.e., include FR-I RGs and BL Lac objects the scatter in- 
creases due to peculiarities of these objects. The scatter further 
increases if we included objects FKM classify as high state ob- 
jects. Thus, the reduced scatter supports the classification of 
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Table 3. Effect of the observing flux limits on correlated data: Most radio fluxes in our sample are obtained with the VLA, which 
can detect 0.1 mJy within a 10 minute snapshot. Chandra can go as deep as 10~ 14 erg/s/cm 2 in a 1000 s observation. 



Radio limit X-ay limit £r £m b x 



1.4 ±0.11 -0.86 + 0.15 -4.84 + 3.2 

0.5 10~ 13 1.37 ±0.11 -0.80 ±0.15 -4.21 ± 3.2 

5 10~ 13 1.34 ±0.10 -0.77 ±0.13 -3.67 ±3.1 

0.5 10- 12 1.36 ±0.12 -0.76 ±0.16 -3.94 ± 3.4 

5 10 12 1.27 ±0.10 -0.64 ±0.12 -1.75 ±2.8 



AGN classes by FKM and suggests that there is a difference 
between LH state AGN and HS state objects. 

3.3. Interpretation in the context of the proposed 
models 

In Fig. |4] we show the % z maps of the different samples and 
the predictions of the jet model and one disk/jet model. For 
the 'jet only' model the predicted values are: £r = 1.38 and 
£ M = -0.81 (FKM) while the values for the 'ADAF/jet' model 
are £ R = 1.64 and <f M = -1.3 (MHDM). The exact value for 
a disk/jet model depends on the solution used for the accretion 
flow. 

The full KFC sample contains some sources, for which we 
used optical fluxes to derive our equivalent X-ray luminosi- 
ties. Thus, the 'ADAF/jet' does not have to be valid. If we ig- 
nore this, the full KFC sample seems to favour the 'jet only' 
model as, according to the confidence region, we can rule out 
the 'ADAF/jet' possibility by more than 3 cr. However, the ex- 
act fit values depend on the choice of the sample and the used 
assumptions for the statistical model. 

In contrast to the claims of iHeinzl J2004I) . the correlation 
is in agreement with a simple jet model in the regime where 
radiative cooling is not important. The choice of sources mini- 
mizes the effect of radiative cooling and the synchrotron cut-off 
as discussed in Sec. 12.51 but it can not be ruled out that cooling 
has some effect on this sample. 

The correlation with the least scatter is found for the sub- 
sample of the KFC sample containing only LH state XRBs, Sgr 
A* and LLAGN. Its confidence region is shown in the middle 
row of the figure. As this subsample consists of 'real' radio/X- 
ray data, both models claim to be valid. However, the fit is not 
in good agreement with both models. The 'jet only' model is 
disfavoured with « 3<x while the 'ADAF/jet' model is even 
stronger rejected. The remarkable low scatter in this subsample 
<Ti„ r = 0.10 supports the idea that LH state XRBs and LLAGN 
are indeed associated. 

In the light of the jet model, it may be that the deviation 
of the fit from the predicted value is due to synchrotron cool- 
ing and the synchrotron cut-off. Even though we have designed 
this sample to minimize their effect, it may still play a role in 
some of the objects. To analyze this effect, one would have to 
compare the X-ray spectrum with the assumed hard power-law 
and take the spectral index into account. If one finds signifi- 
cant deviations in the spectral index one will have to resort to a 
more complicated study where the measured spectral index or 



a more complicated model of the SED is taken into account in 
the fitting. Another possible explanation for this deviation, is 
that we do not treat the coupled errors of the XRB data points 
correctly. As we treat all errors as independent we do not con- 
strain the correlation index £r from the XRBs as well as we 
could. If we set this index to £r = 1 .4 and only fit £m we find 
£m = -0.85 for the KFC subsample containing LLAGN and 
XRBs. See also Fig. [3] This is roughly in agreement with the 
value predicted by the jet/synchrotron model of £m = -0.81. 

For a 'disk/jet' model the discrepancy may be due to the 
disk model used. Besides the discussed ADAF solution, one 
can use any other accretion flow model to create the X-rays. 
This can change the prediction considerably, see e.g., MHDM. 
Thus, the discrepancy of the fit compared to the model predic- 
tions can not rule out any of the suggested models with cer- 
tainty, but - in contrast to earlier claims - it does not support 
them either. 

To compare the edited MHDM sample with the models, we 
now have to note that the 'jet only' model claims to become in- 
valid for the X-ray emission of high state objects like Quasars. 
Interestingly, the edited MHDM sample including high-state 
objects indeed disfavours the 'jet only' model and is in agree- 
ment with the 'ADAF/jet' model (1.2 cr). However, also for this 
sample the selection effects are dominating the exact fit value 
as well. 

3.4. The conspiracy 

We have seen that the fundamental plane in the two described 
incarnations is only slightly different. For both samples we find 
only slightly different parameters, and the scatter seems to be 
less in the case of the KFC sample. The MHDM sample con- 
tains low luminosity objects as well as bright quasars. However, 
there do not seem to be obvious outliers. 

The radio emission is usually attributed to the jet. For the 
higher observation frequencies there are objects in the two sam- 
ples that are clearly jet dominated and others for which the ac- 
cretion flow will be the dominant part at higher frequencies. 
The clearest examples for synchrotron emission are the BL Lac 
objects in the KF C sample but, also fo r FR-I RG this origin is 
well established JChiaberge et alll9 99). On the other hand, the 
X-ray emission from radio quiet Quasars is very likely not syn- 
chrotron emission. Nevertheless, even though the inclusion of 
the Quasar subsample changes the correlation and increases the 
scatter, they do not drop off the correlation like HS state XRBs. 
Likewise, if the X-rays of LH state XRBs and LLAGN are ere- 
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ated by the accretion flow, why do the BL Lac objects and FR- 
I RGs still follow the fundamental plane? There seem to be 
a "fundamental plane conspiracy": even though the emission 
processes are different the objects all lie near the fundamental 
plane. 

Within the jet model one can explain part of the conspir- 
acy by the different emission processes (see Fig.|2}. If we ob- 
serve a source at a frequency after the synchrotron cut-off, the 
Compton branch takes over. The inverse Compton emission can 
in many sources reach the values one would find if one extrap- 
olates the synchrotron power-law to X-ray frequencies. It will 
mainly increase the scatter in the correlation. 

4. Conclusions 

In the previous sections we have reconfirmed the existence of 
the fundamental plane of accreting black hole in the black hole 
mass, radio and X-ray luminosity space. We find that the re- 
sult of a statistical analysis of the radio/X-ray correlation de- 
pends strongly on the assumptions of the distribution and mag- 
nitude of the measurement errors and the intrinsic scatter. The 
measurement uncertainties have been taken from the literature. 
The unknown intrinsic scatter, e.g., the scatter due to relativis- 
tic beaming, non-simultaneous observations or source peculiar- 
ities, has been parameterized and estimated for the observed 
samples. 

Using this refined method we compared the proposed 
radio/X-ray correlations of MHDM and the improved KFC 
sample based on FKM. Both samples differ in their source se- 
lection: while the KFC sample tries to include only sources be- 
longing to the low/hard state, the MHDM sample includes all 
kinds of AGN. Also the observing frequencies differ for some 
sources, as the KFC sample uses extrapolated optical fluxes for 
FR-I RGs and BL Lac objects. 

The best fit values of both samples depend on the relative 
number of sources in each class of objects, e.g., the relative 
number of quasars or FR-I RGs compared to LLAGN. This 
can be understood if different physical processes are dominant 
in the different classes, e.g., the emission from LLAGN may be 
due to the jet, while for quasars the disk emission dominates. 
The confidence regions do not reflect this problem and have to 
be viewed as a lower bound on the errors of the parameters. 

The best fit values found for the KFC sample are £r = 
1.41+0.11 and^M = -0.87 +0.1 4 while we find for the MHDM 
sample^ = 1. 74 + 0.20 and £ M = -1.35 +0.27. Thus, the KFC 
sample suggests a simple uncooled 'jet only' model while the 
MHDM sample favours the 'ADAF/jet' model. However, the 
selection effects are very hard to control. 

The KFC sample seems to be a more homogeneous sam- 
ple, as it has a lower intrinsic scatter. The fundamental plane 
for the subsample containing only LLAGNs and XRBs is sur- 
prisingly tight with a scatter of cr ml = 0.12 dex, while the full 
sample has <j m , = 0.38. Compared to this, the MHDM sample 
has a higher intrinsic scatter of cr mt =s 0.6 dex. This supports 
the AGN classification of FKM in low/hard and high/soft state 
objects. 

In general, the fundamental plane of black hole activity is 
confirmed by our analysis. With a careful control of a homoge- 



neous source selection (high-state versus low-state), the scatter 
can reach rather low values. This promises a wider application 
of the "fundamental pla ne" in other contexts (see e.g.. lMerlonil 
l2004HMaccaronel2 005) and calls for improved radio and X-ray 
surveys in the future. 
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